Luttinger liquid properties of the steady state after a quantum quench 
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We study the dynamics resulting out of an abrupt change of the two-particle interaction in models 
of closed one-dimensional Fermi systems: the field theoretical Tomonaga-Luttinger model and a 
microscopic lattice model. Using a nonperturbative approach which is controlled for small two- 
particle interactions we are able to reach large times allowing us to access the properties of the 
steady state of the lattice model. Comparing those to the exact solution of the full dynamics in the 
Tomonaga-Luttinger model we provide evidence for universal Luttinger liquid behavior. We show 
that a single impurity leads to open and perfect chain fixed points in the steady state similar to the 
ones found in equilibrium [Phys. Rev. Lett. 68, 1220 (1992)]. 
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With the rapid progress in the preparation and mea- 
surement techniques for isolated cold gases [l| investi- 
gating the fundamental questions of if and how a closed 
quantum many-body system prepared in a nonequilib- 
rium initial state approaches a stationary one is within 
experimental reach. Studying the physics of the steady 
state itself is of particular interest if it is nonthermal [^ , 
that is expectation values of observables differ from those 
computed using a canonical statistical operator with the 
temperature fixed by the excess energy. Onc-dimcnsional 
(Id) interacting Fermi systems are promising candidates 
for realizing such unusual nonequilibrium states as in 
many of those the dynamics is restricted not only by en- 
ergy conservation but by several additional conservation 
laws [3|. An often studied protocol is an abrupt param- 
eter quench. The system is prepared in the canonical 
ensemble of an initial Hamiltonian H(Ui). The time evo- 
lution is performed with H{U{), C/f ^ Ui. After taking 
the thermodynamic limit a finite subsystem [J] or local 
observables [5| might become stationary and describable 
by a time independent statistical operator. In our case 
the quenching parameter U is the amplitude of the two- 
particle interaction. 

Interacting Id quantum systems are exceptional for a 
different reason as well [6^, Jl . The equilibrium low-energy 
properties of generic models are not described by a Fermi 
liquid theory. Instead of being fermionic quasi-particles, 
the gapless elementary excitations are of bosonic nature. 
Systems of this type fall into the Luttinger liquid (LL) 
universality class [8[ . The observation that the field the- 
oretical Tomonaga-Luttinger (TL) model is the infrared 
fixed point of a large class of models under a renormaliza- 
tion group (RG) flow lies at the heart of the equilibrium 
LL concept [9| . For the case of spinless models the fixed 
point is characterized by the velocity v of the collective 
excitations and the LL parameter K entering exponents 
of algebraically decaying correlation functions. Both de- 
pend on the parameters of a given model, in particular 
on U . Alternatively, the two-particle interaction might 



open a gap, a situation we do not consider here. 

One of the hallmarks of LLs is their sensitivity towards 
inhomogeneities. For repulsive interactions with K < 1 
the density response function of a LL diverges as [10| 
x((7, oj = 0) ~' |(7 — 2A:fP^~^, with the Fermi momen- 
tum fcp. This indicates that even a single weak impu- 
rity acts as a strong perturbation. In fact, the homo- 
geneousperfect chain fixed point turned out to be un- 
stable llMl3l|. The system flows towards an open chain 
one with dramatic consequences for the equilibrium low- 
energy properties; e.g. the linear conductance vanishes 
as G - T'^K-^-2 fQj. temperature T ^ 0. 

In this Letter we achieve two distinct goals, (i) We 

show that the steady state of a microscopic Id lattice 

model after an interaction quench is characterized by the 

same power laws as the TL model after a similar quench 

14l4l9{ with the K taken for the respective model. As the 



concept of RG irrelevance cannot directly be transferred 
to nonequilibrium [16, 19, 2Qj this is a highly nontrivial 
result which complements earlier indications of LL uni- 
versality in the time evolution towards the steady state 
[2l|. (ii) We study the physics of a single impurity in 
the nonequilibrium steady-state LL. For K < 1 a sta- 
ble open and an unstable perfect chain fixed point is 
found, in analogy to equilibrium ll[. We determine the 



power-law scaling of the conductance as a function of 
T. The data for different impurity amplitudes fall on a 
one-parameter scaling curve indicating the absence of any 
additional fixed points. The dynamics of the TL model is 
solved exactly using bosonization [y, [Zl, [2^ . To study the 
time evolution of the lattice model we use an approximate 
functional RG [23| based approach which was originally 



developed for open quantum systems [2J]. We establish 
it as a tool for studies of closed ones as well. This il- 
lustrates the remarkable flexibility of the functional RG 
approach to quantum many-body physics. For small two- 
particle interactions this essentially analytical approach 
allows controlled access to time scales large enough such 
that the physics is dominated by the steady state. This 



complements numerical approaches which provide 'exact' 
results at small t but abru ptly become unreliable beyond 
a characteristic time scale |25l428| which might be smaller 
than the one on which the steady state is reached 21 1. 

Bosonization — The starting point of our investigation 
of LL universality is the exact computation of the desired 
observables and correlation functions within the (spin- 
less) TL model using bosonization [6|, [Zl, l22|. Starting 
out from the Id electron gas the TL model is obtained by 
linearizing the single-particle dispersion and keeping only 
the marginal two-particle scattering terms. In contrast to 
earlier studies on interaction quenches in the TL model 
lj-|2l|, l29l - l31| we consider the one with open boundaries 
at a; = and x = L [ll|, [S^l ■ This allows us to distin- 
guish between bulk and boundary LL exponents 33ll34|. 
While the former are quadratic in the two-particle inter- 
action the latter are linear. The model is given by the 
Hamiltonian 



OO r _, 



(1) 



with the Fermi velocity v-p, kn = nn/L, the two-particle 
potential u(fc), and bosonic operators 5,/ associated 
with the density of the fermions. Employing a Bogoli- 
ubov transformation -ffxL can be written as a diagonal 
quadratic form in cigcnmodes with (bosonic) ladder op- 
erators a„ and energy cjn = VFkn\/l + u{kn) / {ttvf) ■ 
The mode occupancies al^a„ are independent integrals 
of motion. To keep the formulas transparent we here 
take the nonintcracting canonical statistical operator p'^ 
as the initial state; in [35| we describe the changes when 
starting in the canonical state with Wi(fc) > 0. 

Using the Bogoliubov transformation and bosonization 
of the field operator [gI, [ZI, 122, |32| it is straight forward to 
derive closed analytical expressions for the time depen- 
dence of a variety of observables and correlation func- 
tions. We exclusively consider quantities which become 
stationary and aim at their steady state values. Those 
are obtained by computing Tr [p°e*^*Oe~*^*] and af- 
ter performing the thermodynamic limit taking the limit 
t ^ OO, where O stands either for an observable or the 
operator product defining a correlation function. We 
verified that the same expectation values follow from 
taking limi^oo Tr [pggeO], with the statistical operator 
Pgge = exp{-X]„>o^n"l"rJ/-^GGELof a iionthermal 
generalized Gibbs ensemble (GGE) [3l. Il4l lla. llTl Il9l |36| 
and ^GGE determined such that TrpcGE = 1- The rjn 
are fixed by requiring Tr [a\a^pQOF\ = Tr [a]^Q;„p°] . 

Despite being described by the density matrix pgge 
the steady state of the TL model is 'critical', that is char- 
acterized by LL power laws. This is known from quen ches 

illii 



obscrvablc/corrcl. funct. variable cq. cxp. steady-state exp. 

access density An x —K ^{^ + l)/2 

local DOS p u) R-^ - 1 (R-^ - l)/2 

bulk x at tj = q - 2fcF 2{K - 1) K'^ - 1 



TABLE I. Equilibrium and steady-state scaling exponents. 



can all be expressed in terms of the models LL param- 
eter K = [1 + u{0)/{TrvF)]~^^^ but differ from the cor- 
responding known equilibrium exponents (third colum) 
jd-lllj. In the next section we directly compare the de- 
cay of the densities Friedel oscillations off the boundary 
in the steady state of a microscopic lattice model with 
the TL model prediction (first row of Table |T| . 

The scaling behavior of the bulk density response 
x{q, oj = 0) (last row of Table|l| allows us to make predic- 
tions for how the steady state reacts to a single impurity. 
For repulsive interactions K^ — 1 < and x diverges for 
q — ^ 2fcF. As in equilibrium even a weak single local im- 
purity strongly disturbs the homogenous system. When 
applying an infinitesimal bias voltage V across the impu- 
rity the conductance G = dP^ /dV (with the stationary 
current P^) is expected to scale as Go - G{T) -^ T^ ~^, 
with the constant homogenous chain conductance Gq. In 
the following we denote G as the steady-state linear con- 
ductance. The noncquilibrium nature of the considered 
state stems dominantly from the interaction quench and 
not the small \V\. The power law holds as long as the 
right hand side stays small, that is for not too small T. 
Using the language of equilibrium RG this indicates that 
the perfect chain fixed point is unstable. In contrast, the 
open chain one is stable as follows from the scaling of 
the local density of states (DOS; second row of Table |T|. 
Fermis Golden Rule-like arguments lead to the vanishing 
of the tunneling conductance G ~ T^ ~^ across a weak 
link connecting two semi-infinite chains with suppressed 
DOS. These arguments do not rule out intermediate im- 
purity fixed points. As in equilibrium |ll| excluding those 
requires a method beyond perturbation theory in the im- 



14 



and 



in the translational invariant TL model 

intimately linked to the mode index n dependence of 

rjn/iOn- Table H] contains the T = scaling exponents in 

the nonequilibrium steady state of interest to us. They in standard second quantized notation. In equilibrium 



purity strength and the weak link hopping |12l . Il3j . Pro- 
vided the concept of LL universality holds for the steady 
state we expect to find those weak and strong impurity 
scaling laws of G{T) for a lattice model with the K of the 
model considered. Using a one-parameter scaling ansatz 
and our nonperturbative RG we investigate the possibil- 
ity of intermediate impurity fixed points. 

Microscopic model — We consider the lattice model of 
spinless fermions with nearest-neighbor hopping J as well 
as interaction U and open boundaries terminating the A^- 
site chain given by 

Hlm{U)=J2 [-^4s+i +H-C. + Uc]c/^^,c^^,] (2) 
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FIG. 1. (Color online) Functional RG data for the time evo- 
lution of rij (i) — u at half filling v = 1/2 after a quench in the 
interaction amplitude to U/J = 0.5 for A'^ = 10^, T = 0, and 
different sites j. Left inset: Friedel oscillations induced by 
the boundary and the propagation of a main signal from the 
boundary to j ~ vt for two values of t. Right inset: Velocity 
of the main signal for different U (symbols). The exact Bethe 
ansatz v (line) is in excellent agreement with our data. 



the model is (A) Bethe ansatz solvable and (B) shows 
universal LL physics with K and v exactly known [8[. 
It is commonly believed that because of (A) the steady 
state after an interaction quench is described by a GGE 
but the corresponding pgge was neither constructed nor 
was a proof of its existence given. Our analysis does not 
rely on any such assumption. When discussing the impu- 
rity physics Hj^M is supplemented by a hopping impurity 



FIG. 2. (Color online) Functional RG data for the effective 
exponent of the Friedel oscillations of the steady-state density 
nf for u = 1/2, A^ = 10^, and T = determined by taking the 
log-derivative. The predictions from the TL model with the 
exact lattice model K (dashed lines) are consistent with our 
results. The inset shows In^' — !/| skv = 0.375 and U j J — 0.25 
(symbols). The line is the TL model prediction. 



the left signal passed through and the right one did not 
enter yet the density becomes stationary. The physics 
for t w 10^/J and j up to 0(10^) is thus dominated by 
the steady state and does barely suffer from finite size 
effects. The two main signals propagate with the LL ve- 



locity V [37 



ffi 



imp 



'^'^N/2^N/2+l 



- 39[ , with our method providing an excellent 
approximation to the exact v (right inset of Fig. [1}. 

In analogy to the interacting ground state (see e.g. [12|) 
the stationary density nf shows Friedel oscillations with 
frequency 2fcF (see the j < 80 region of the lower part 
H.c. of strength h e [0, J] located of the left inset of Fig. [T]). We next analyze their decay. 



in the middle of the chain. To compute the time evolution 
of the density rij as well as the conductance we use an ap- 
proximate functional RG [23| based approach originally 
developed to investigate the nonequlibrium dynamics of 
open systems [2J]. A related method was earlier shown 
to capture the equilibrium LL properties of inhomoge- 
neous lattice models 12|, |l3| including the characteristic 
power-law scaling, with exponents agreeing to the exact 
ones to leading order in U. Our approach is essentially 
analytical but in its last step requires the numerical so- 
lution of a set of 0{N-Nt) coupled differential equations, 
with the number of time steps Nt , and the computation 
of matrix exponentials of iV x A^-matrices [35| . 

In Fig. [1] we show the access density nj{t) — v for fixed 
j starting out of the noninteracting impurity free ground 
state (T = 0) [35[. Despite the open boundaries for fill- 
ing V = 1/2 the initial state has a homogeneous density 
nj{t = 0) = 1/2 3j| . We can reach times of the order of a 
few 10^/ J which has to be contrasted to other commonly 
employed methods which become unreliable for times of 
the order of 10/ J [21j. As shown in the main plot and 
the left inset a prominent signal originating from the left 
boundary travels through the system. A similar one is 
generated at the right one. For a spatial region in which 



Figure [5] shows the log-derivative of \nf — 1/2|, that is 
an effective exponent. The dashed lines is the predic- 
tion from the TL model of Table U with the exact lattice 
model K . Our data are consistent with a power-law de- 
cay and the TL model exponent. This finding is our first 
indication of LL universality of the steady state. The 
differences between the exact exponent and our result 
is of order t/^. On the right hand side of our RG flow 
equations we do not fully capture terms ~ U^ and thus 
control exponents only to order U . The discussed behav- 
ior is not restricted to the case of half filling. The inset of 
Fig. [2]shows \Tvf — v\ for v = 0.375 on a log-log scale and 
the corresponding TL model prediction as the envelope. 
To compute the steady-state linear conductance of the 
lattice model we take the canonical density matrix pj? 
(with T > 0) corresponding to Hi.m{U = 0) + i/imp as 
the initial state. The time evolution is performed with 
HhuiU > 0) + i/imp supplemented by onsite energies 
V/2 i-V/2) for all j < N/2 {j > N/2). Thus not only 
the two-particle interaction but in addition a small bias 
voltage V is abruptly turned on at t = 0. The current /(t) 
across the impurity bond is computed [35[ . Following the 
same reasoning as for the density / becomes stationary 
for t of the order of 10^/ J. We take V to be the smallest 
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FIG. 3. (Color online) Functional RG data for the tempera- 
ture derivative of the steady-state linear conductance of our 
lattice model for a weak impurity h/ J = 0.05, filling u = 1/2, 
and A'^ = 10'^ (symbols). The prediction from the TL model 
with the exact lattice model K is shown as the solid lines. 
To emphasize the sizable differences in the exponents we have 
added a power law with the U / J = 0.1 exponent as the dashed 
line to the U / J = 0.5 data. 



energy scale of the problem (typically V = 10 ^J) to 
ensure that we are in the linear regime /^* = GV. For 
T > J we find independent of h, U, and u, G{T) - T'^ 
[l3| which is a band effect (see the inset of Fig. U). The 
universal scaling of the conductance discussed in the last 
section can only be expected for T <^ J . 

We first analyze the case of weak impurities. To elimi- 
nate the constant Go we take the derivative of G with re- 
spect to T. Based on our above considerations we expect 
to find dG/dT - T^'-^; see Fig. O Over roughly one 
order of magnitude the functional RG data follow the TL 
model prediction with the exact K of our lattice model. 
The deviations for T/J > 0.1 indicate the crossover to 
the nonuniversal G{T) ^ T^^ regime. The deviations 
for T/J < 0.01 have two reasons. As discussed above 
the scaling only holds as long as T does not become too 
small. Furthermore, the energy level spacing Sn ~ vp/N 
(= 2 • 10"'^ J for the parameters of the plot) is an energy 
scale of the problem which cuts off any universal behavior 
12, ll3| • This is an artefact of our treatment of finite sys- 
tems. For small /i and T/J e [0.1, 0.005], Go-G'(T) < 1. 
Our analysis thus requires very accurate data. To min- 
imize the error due to small residual oscillations of I{t) 
present even for t of the order of 10^/ J we averaged the 
data at large t over a small time interval. 

In the inset of Fig. |4] wc present our results for G{T) 
across a strong impurity {h/J = 0.9). Even without any t 
averaging our data are accurate enough such that the log- 
derivative, that is the effective exponent, gives a smooth 
curve. The data clearly show the crossover from the 
nonuniversal T~^ behavior at large T to the TL model 
prediction T ^ at low ones down to the cutoff scale 
5n- In the limits of strong and weak impurities our re- 
sults for the linear conductance of the lattice model thus 
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FIG. 4. (Color online) Functional RG data for the one- 
parameter scaling of G{T) of the lattice model for U/J = 0.5, 
V — 1/2, and A'' — 10'^. Different symbols stand for differ- 
ent h increasing from left to right. No fixed point in between 
the perfect and the open chain ones exist. The crosses were 
computed for U/J — 0.85, v = 1/4 giving the same K. The 
scaling function thus depends on U and u only via K [ill ]. 
The prediction of the TL model for a:: — )■ and a:: — > oo with 
the exact lattice model K is shown as dashed lines. Inset: the 
effective exponent of G{T) for a strong impurity h/J = 0.9 
determined by a log-derivative. Dashed lines show the pre- 
diction of the TL model with the exact lattice model K. 



agree to the TL model prediction providing the second 
indication of LL universality of the steady state. 

We finally show that in the steady state of the lat- 
tice model no fixed point in between the perfect and the 
open chain ones exist. To this end we compute G{T) for 
a variety of h at fixed U and v. Using a one-parameter 
scaling ansatz of the form G{T) = G{x) [ll|, ll3|, with 
X — {T/sY^ -i)/2 and the nonuniversal scale s{U,i',h), 
all data can be collapsed on a single curve continuously 
connecting the weak (x — >■ 0) and strong {x -^ cx)) impu- 
rity fixed points; see the main plot of Fig. [H We conclude 
that the impurity physics of the LL steady state is similar 
to the equilibrium one |ll| . 

Outlook — We provided evidence that the steady state 
of an interacting Id lattice model after a quantum quench 
shows LL universality and that the impurity physics of 
the LL steady state is similar to the one in equilibrium. 
Both should be confirmed by studying other models of 
interacting fermions using alternative methods. In par- 
ticular, it would be interesting to investigate if the 'criti- 
cal' LL physics can be found in the steady state of mod- 
els which are not Bethe ansatz solvable ('nonintegrable'). 
For the asymptotic time dependence indications of uni- 
versal LL power laws were found 2l|. Studying 'nonin- 
tegrable' Hamitonians obtained from the one Eq. ([2]) by 
adding further terms is beyond the current possibilities 
of our method as wc would have to give up the crucial 
requirement of small nearest neighbor interactions. 
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BOSONIZATION 

To compute observables and correlation functions in 
the steady state of the Tomonaga-Luttinger (TL) model 
with open boundaries after an interaction quench we use 
'open boundary bosonization' for the Hamiltonian and 
the field operator ^\^. We are interested in the scaling 
behavior with all energy scales send to zero and all length 
scales send to infinity. Thus subtleties resulting out of the 
momentum dependence of the two-particle potential u{k) 
become irrelevant [SJ, |6| and the ultraviolet regularization 
can be implemented at will. To illustrate the procedure 
we consider the density n(x, t). We first study the quench 
from Ui{k) = to Mf(fc) = u{k). For simplicity we focus 
on temperature T = 0. The density is given by the Green 
function 



iGt{x,x) = (vac(&)|e*^^^V^(a;)V'(a;)f 



-iHTLt 



vac(5)) 



with the noninteracting ground state |vac(6)) which cor- 
responds to the vacuum with respect to the &'s (in the 
same notation as in the main text). The fields ?/'^^'(x) 
contain the open boundary conditions. Using auxiliary 
fields V''^-* (x) which arc identical to the ones obtained for 
periodic boundary conditions [7h9[ and are e.g. given in 
Eqs. (18)-(20) of Ref. [y], the Green function reads 



iGt{x,x)-- 



Ir 



;vac(6)| e'"-^'^*i>'i{x)^{x)e-'"-^'^* |vac(6)) 

;vac(6)| eJ^-^'^^ij)'' {-x)iP{-x)e-'"-^'^* |vac(6)) 
;vac(6)| e''-^TL*V'^(^)^(-a;)e-*-^^L* |vac(&)) 



{Yac{b)\e'"-^''*iP^-x)ip{x)e 



-iHtl* 



vac (6)) 



Those expectation values can be computed following the 
usual steps [7H9[ which involve the multiple use of the 
Bogoliubov transformation 6„ = c{kn)ct^ + s(fc„)ajj and 
the Baker-Campbell-Hausdorff relation. The coefficients 
c{kn) and s(fc„) depend on the two-particle potential 
u{kn) and are e.g. given in Eq. (9) of Ref. [y]. The first 
two terms of the Green function provide the homoge- 
nous background density while the latter two oscillate in 
space with frequency 2fcF — they contain the Friedel os- 
cillations induced by the boundaries. After taking the 
thermodynamic limit and the limit i — > oo we find for 
the steady-state access density 



At zero momentum the coefficients of the Bogoliubov 
transformation can be expressed in terms of the mod- 
els Luttinger liquid (LL) parameter K given in the main 
text as 



s^{0) = \{K + K-' 



2), c^{0)^^{K + K-'+2).{2) 



Using those relations the exponent of Eq. ([T]) can be writ- 
ten as -(A'2 + f )/2 (sec Table I of the main text). The 
other scaling exponents of the last column of Table I of 
the main text can be obtained in a similar fashion. 
In any (equilibrium) LL the LL parameter is given by 

0{U^), (3) 



K ^l-U 



where C/ is a dimensionless measure for the interaction 
strength; e.g. U = u{Q) / {2Trv-p) for the TL model and 
U ~ U[l — cos(2A:F)]/[27rJsin(A:F)] for our lattice model 
at filling v = kp/n. Using this expansion it is evident 
that all scaling exponents discussed by us (sec Table I of 
the main text) have a leading order contribution in U. 
This is crucial as within our approximate treatment of 
the lattice model we control exponents only to leading 
order 0. 

We next briefly discuss the case when starting in the 
ground state of H^h with Ui{k) > and performing the 
time evolution with a different interaction Uf{k) > 0; all 
quantities depending on the interaction strength acquire 
indices i or f. We have to consider two Bogoliubov trans- 
formations [corresponding to the transformations from 
zero interaction to Ui/{(k)] and two sets of eigcnmode 
ladder operators. The initial state is the vacuum with 
respect to one of those while the Hamiltonian with which 
the time evolution is performed is a diagonal quadratic 
form in the other one. Repeatedly applying the Bogoli- 
ubov transformations and the Baker-Campbell-Hausdorff 
relation gives for the steady-state access density 



An'^*(a;)-5 
with the scaling exponent 



cos(2A:Fa;), 



(4) 



7 = K{ 



1 



1 + -(/M + Kr' - 2)iK{ + K-' H 

o 

\{K,~Kr^){K,-K-') 



2) 



(5) 



l^Tf^{x) - X' 



[c\0)+s\0)][c{0)+sm 



cos{2kpx). (I) 



For Ki = 0, that is if we start in the noninteracting 
ground state, it becomes equal to [Kf 4- 1)/2 (see Table I 



of the main text). Using the expansion Eq. (jS]) it is easy 
to see that 7 and {K'^ + l)/2 agree to leading order in 
the two-particle interaction, that is Ui only contributes 
to order Uf and C/fC/j or higher. To leading order in the 
interaction strength the scaling exponent of the density 
is thus exclusively given by C/f. The same holds for the 
other exponents considered by us (see Table I of the main 
text). In our computations for the lattice model we con- 
trol exponents only to leading order which explains why 
in the main text we exclusively consider quenches out of 
the noninteracting ground state. 



FUNCTIONAL RENORMALIZATION GROUP 

In the main text we use the functional renormalization 
group (RG) approach [lO[ to study the relaxation dynam- 
ics and the steady state of a closed many-body system 
described by the lattice model of spinlcss fcrmions with 
nearest neighbor hopping and interaction. Compared to 
the functional RG approach to time evolution for open 
quantum systems introduced in Ref. 



with e being the N x N matrix with entries e^ and n^. 
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some minor 

amendments need to be made. Those arc outlined next. 
We can treat Hamiltonians of the form 



H — Ho + Hint, 


(6) 


ij 


(7) 


-'- \ ' - ft 
int = 7 7 ^UijklC^C^ClCk 


(8) 



ijkl 



written in standard second quantization. Here uij^i is 
the antisymmetrized two-particle interaction. The in- 
dices i,j,... stand for the quantum numbers, e.g. the 
N Wannier states in our lattice model. The Hamiltonian 
Eq. (2) of the main text falls into this class of general 
Hamiltonians. We assume an initial density matrix 



al^ij I 



^^y■\(.T..,P^34<^3 



(9) 



which allows for the application of Wick's theorem. 
In the main text we always choose the noninter- 
acting canonical statistical operator p^ = p^ = 
exp(— /37Jo)/ Tr[exp(— /3_ffo)]- The two-particle interac- 
tion is treated approximately using functional RG in its 
lowest order truncation llOl . We introduce a cutoff in the 
noninteractin g K cldysh [12| Green functions g (as moti- 
vated in Ref. [ll|) by 



^rct,A 



{t,t') = -ie(i - t')e-*<*-*')e-'^(*-*') 



^adv,A 



{t',t)]\ 



f^\t,t') 



^W 



'■'=*'^(i,0)(l-2n)<?^^^'^(0,t') 



(10) 

(11) 
(12) 



Tr 



P°dld^ 



The self-energy is obtained by solving a set 
of coupled differential flow equations 

9As:::;f(t',i)-9As::^;^i',t) 

= - E ^5^*' *) i~'^n^.^[^i. (0) S{t' - t) , (13) 



c^aS^'^ = 0, 




(14) 


with the initial conditions at A = 


: OO 




s:-f^^=°°(t',t) = i^(t- 


-<')]> ]ui,j,j, 


(15) 


S^f=°°(i',i)=0. 


J 


(16) 



The right hand sides of the flow equations contain 



(17) 



with the full cutoff dependent Keldysh component of the 
Green function 

G^'^it, t') = -iG'°^^'^{t, 0)(1 - 2n)G^^^'^(0, t'). (18) 

The star differential operator d^ acts only on the 
free Green function ^rot/adv,A^ ^^^^ ^^^ j^A^ ^^ ^-^^ 
Dyson series expansion 0'^°*/^'^^,^ = c,ret/adv,A _,_ 

^rot/adv,Aj.A^ret/adv,A _^ _ _ ^^^^ ^^ calculate G^-et/adv^A^ 

An approximation to the self-energy of the cutoff-free 
problem is obtained at A = 0. 

We can calculate G^'^{t,t') very efficiently by using 
an iterative procedure. First, we discretize time in steps 
such that during one small step At the time dependent 
self-energy can be set constant. For our results shown 
in the main text we made sure that Ai is always chosen 
small enough such that further reducing it does not lead 
to any changes visible on the scale of the respective plots. 
We use 



G"'''^it,t')G"''^^{t',t") 

= -ie(t - t')e{t' - t")G"'*''^(i, t") 

it 



(19) 
G"'''^it,t') ^ [G"'^"^^(i',i)]\ (20) 

to write G'^°^{t,t') as a product of Green functions 



G''°^'^{ti + At,ii) = -■-■p-^[<^+^"' (*i)l^*^-*A(t-*') 



-le 



(21) 



where fi''''*''^(ii) is the self-energy time averaged over the 
interval (ti,ii -I- Ai). The interacting Keldysh Green 
function G^'^{t, t) can then be found iteratively employ- 
ing 

G^'^{t + At,t + At) 

= G''<=*'^(t + At, t)G^'^(t, i)G''<=*^'^(t, t + At), (22) 
GK■^(0,0) = -^(l-27l). (23) 



In every time step two matrix exponentials oi N x N 
matrices have to be performed and multiplied with the 
Keldysh component of the Green function of the previous 
one. This renders the problem a natural candidate for 
graphics processing unit (GPU) supported algorithms. 
We use such to compute the results shown in the main 
text. The number Nt of time steps needed to obtain sufB- 
cient accuracy (and resolution) as well as to reach times 
which are large enough such that the physics is domi- 
nated by the steady state enters the number of equations 
to be solve. For the Hamiltonian considered in the main 
text we solve sets of {3N — 2)Nt coupled differential equa- 
tions. Due to the nearest neighbor structure of the in- 
teraction 3A^ — 2 components of the self energy flow for 
each of the Nt time steps. Typical numbers considered 
are N = 1000 lattice sites and Nt = 1600 time steps. 



A''. Furthermore, we use 

lit) = -I {[H,NL]it))p„ = {h- J)G</2j^/2+i(t,i) + c.c. 

(26) 
with G^,^j^,2,^{t,t) being the equal-time lesser Green 
function of the interacting system [12|. The functional 
RG method used here provides an approximation for this 
given by 



G<''^=°(M) 



g:^'^^°(m) 



(27) 



Therefore, plugging this lesser Green function into 
Eq. (|26p allows to compute an approximation to the cur- 
rent and from this the conductance. 



COMPUTING OBSERVABLES 

An approximation for the occupancy of site j can di- 
rectly be obtained from the Keldysh Green function at 
the end of the RG flow as 



n,it)^\-'-Gf^^-\t,t). 



(24) 



To calculate the current flowing from the left to the 
right half of the lattice in our microscopic model of in- 
teracting spinless fermions, we need to determine 



Ht)^-J^{NL{t))pO, 



(25) 



where Nl is defined as NL{t) = J2j=i "j(Oi '^itli '"-j = 
ricj being the occupancy operator of site j [l3| and 
(. . .) denotes the expectation value with respect to p*^. 
For simplicity we assume an even number of lattice sites 
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